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Abstract 

We perform an extensive investigation of the sensitivity to non-vanishing mass in a large 
water Cerenkov detector, developing an analysis method for neutrino events originated by a 
supernova explosion. This approach, based on directional considerations, provides informations 
almost undepending on the supernova model. We analyze several theoretical models from numerical 
simulations and phenomenological models based on SN1987A data, and determine optimal values of 
the analysis parameters so as to reach the highest sensitivity to a non-vanishing mass. The minimal 
detectable mass is generally just above the cosmologically interesting range, m ^ 100 eV, in the case 
of a supernova explosion near the galactic center. For the case that no positive signal is obtained, 
observation of a neutrino burst with Super-Kamiokande will anyhow lower the present upper bound 
on Vt mass to few hundred eV. 


1 Introduction. 

Upper bounds on the electron (anti-) neutrino mass at the level of few eV are obtained from /3 decay experiments 
l]^, and have been confirmed by data on neutrinos from S'iV1987A, collected by Kamiokande-II [Q and IMB ||] 
detectors, see e.g. IJ- Concerning the masses of and i^t, the present experimental bounds are much less 
compelling: ^ 170 keV, m ^ 24 MeU ||p. Therefore, any approach capable of significantly lowering these 

limits is highly desirable. Particularly, an observation sensitive to neutrino masses in the range 10 — 100 eV 
would be extremely important for cosmological implications, see |^. Previous analyses |Q, |j claimed that this 
goal might be reached with new generation larger detectors, for a supernova near the galactic center. 

The questions we address in this paper are the following: 

1. Which assumptions about neutrino emission are critical for observing the effects of a non-vanishing mass? 

2. How the analysis procedure can be optimized in order to reach the highest sensitivity to neutrino masses? 

3. Which masses for and Vr can be actually explored with a detector such as Super-Kamiokande? 

In general, the determination of zz-masses by means of a supernova observation is based on the comparison 
between the detected signal and the expected one for the case m = 0. A non-vanishing mass, in fact, would 
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reflect in an overall delay of the events number distribution, and also in a spectral distortion of the signal, due 
to the energy dependence of the flight-time for massive neutrinos: 


At (e^) 


D / TOC^\^ 
2c V Cl/ / 


( 1 . 1 ) 


where D is the supernova distance, c is the light speed, and tv is the neutrino energy. The greatest difficult in 
applying this procedure is in the dependence on the supernova model assumed to calculate the expected signal. 
Moreover, when interested in and Vt (hereafter we will refer to these neutrinos, and their antiparticles, as 
Vi), there is an additional problem: it is impossible to isolate their signal from that of Ve and Ve, because all Vi 
detection processes (namely neutral current reactions) are allowed for electron neutrinos too. 

On the other hand, charged current interactions (only possible for v^ and v^ at the energies of interest) might 
be used to clearly identify part of the electron neutrino signal. One can thus conceive the following strategy: 


1. split, by use of a distinctive signature, the overall set of events in two classes: the first one containing only 
Ve and Vf., the second one involving all neutrino species; 

2. use the first class of events to derive the expected distribution for Ve and i>e- In view of the strict limit on 

we can treat these neutrinos as massless for our purpose (e.g. from rriv^ ^ 15 el^ and equation 1.1, 
we get At ^ 0(1 s) even for the lowest detectable energies, tv — 5 MeV, for a SN near the galactic center); 

3. use this distribution to infer that of v^ and Vr, again for massless neutrinos; 

4. compare the observed signal with that theoretically built in the previous step and look for mass effects. 
We assume that only one neutrino (say v^) is massive and that flavor oscillations do not occur (see ^ and 
references therein for search of Vr mass by use of MSW effects). The comparison should be performed in 
the second class of events {ve + Vi): in such a way we deal with a subset which still includes the whole Vi 
signal, but only a part of the Ve “background”. 


The residual dependence on the supernova model is essentially in step 3: how to connect the Ve and v^ emission 
features with those of Vi. In this paper we extensively investigate this point, by analyzing several theoretical and 
phenomenological SN models and their implications for the extraction of a non-vanishing Vr mass. 

We also discuss how the analysis procedure can be optimized, by a proper choice of the parameters, such as 
energy and time windows where to look for mass effects. 

Concerning the sensitivity to neutrino mass, two points are interesting: the minimal detectable Vr mass and 
the upper limit which can be established on m when no positive signal comes out. Correspondingly, we shall 
address the following questions: 

3a. What is the minimal m value expected to give a signal distinguishable from the statistical fluctuations of 
the massless case? 


3b. What is the minimal m value for which the signal expected for massless neutrinos can be discriminated 
from the statistical fluctuations of the massive case? 


The paper is organized as follows: we develop a method of data analysis, based on directional considerations 
(Section 2) which we use to determine the minimal detectable v^ mass for a set of theoretical supernova models 
(Section 3). The same analysis is repeated for a set of phenomenological models derived from SN1987A data 
(Section 4). We then compare our results with those of ref. |]^ (Section 5). Finally, we calculate, both for 
theoretical and phenomenological models, the mass upper bound obtainable through the same method (Section 
6). Detector characteristics, cross sections and expression for event rates are collected in the Appendix. 
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2 The analysis method: isotropic vs. directional signal. 

Krauss et al. 0 proposed in 1992 a beautiful procedure of data analysis for Cerenkov detectors, based on 
directional considerations. In this section we briefly recall the idea and then we specify our approach to this 
method in terms of: i) the kinematical characterization of neutrino events; ii) the assumptions needed to connect 
the emission features of electron and other flavors neutrinos; iii) the quantitative criterion for the evidence of a 
non-vanishing Vt mass. 

2.1 Generalities. 

Let us consider Super-Kamiokande (SK hereafter), a water Cerenkov detector whose main characteristics are 
listed in the Appendix, see also |^. Neutrino events are originated basically in processes of four kinds: 


1. capture on a proton: 

W j- 

Ve+ P —n + e 

(2.1) 

2. scattering off an electron: 

(-) , _ w.z (-) , _ 

Ve +e —> Ue +e 

(2.2) 


Vi e~ —> Vi + e~ 

(2.3) 

3. capture on an oxygen nucleus: 

z/e + -b e- 

(2.4) 


+ 16 O 167 V + e+ 

(2.5) 


4. coherent scattering off an oxygen nucleus: 

Vx + ^ Vx + x = e,n,T (2.6) 


Expressions for the cross sections are collected in the Appendix. The approximate numbers of expected events 
in SK, for a 5'V1987A-like supernova near the galactic center, are shown in Table 1, see also 0- 


I'eP - 

e^n 

~ 5000 

1 

1 

1 

t 

~ 100 

ViC 

Vie 

~ 2 X 25 

VeO - 

■r Ve 0 

~ 50 

VxO - 

-> VxO 

- 200 (T| 


Table 1: Expected events in SK for a supernova at galactic center. 


Among the four processes, scattering off electrons is the only directional one. For capture reactions and 
interactions with nuclei, the observable particles (e^ or 7 ) are emitted (almost) isotropically, whereas in processes 
(2.2-3) the recoil electron’s angular distribution is strongly forward peaked. Therefore one can define an 
appropriate observation angle 'dfw, along the supernova direction, which includes all of the electrons from v — e 
scattering. 

Here is the required signal splitting: outside the forward cone, nearly any event is originated by v^. and 
The only contribution of other flavors neutrinos is due to (2.6), and can easily be excluded by raising the minimal 
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detected energy up to 7 - 8 MeV |^ . So we can regard this signal as untouched by mass effects, and use it to 
derive the expected distribution in the forward cone for the case m = 0. 

The analysis procedure is sketched in Figure 1. 


2.2 Kinematics. 

A specification of the forward observation angle d/m for the analysis. 

When a neutrino with energy Ci, scatters off an electron at rest which gets kinetic energy Ce, the electron 
scattering angle 'd is determined from energy-momentum conservation: 


cos d(e,y, Ce) 


Ce 

e» -I- 2m, 



(2.7) 


As is a decreasing function of Ce and a rising function of Ci/, the maximum value if max is obtained by taking in 
(2.7) the minimal detected energy emm and taking the limit oo. In this way one has: 


— aCOS 


2me 


( 2 . 8 ) 


For the lowest value of Cmin (he. the energy threshold of SK, eth = 5 MeV) this gives dmax — 24°. 

In Figure 2 we plot the angular distribution of the Vr scattered electrons along the supernova direction, for 
two significative values of the lower cut: Cmm = 5 MeV (SK threshold) and Cmin = 15 MeV (the highest value 
useful for the analysis, as we shall see). Although for these calculations we used the “reference” model to be 
described in Section 3, very similar distributions are obtained from other supernova models. In the same figure 
a plot of 'dmax is given as a function of the lower energy cut. 

Besides this purely kinematic dispersion, a systematic angular uncertainty is also present, due to the detector. 
This uncertainty is a decreasing function of the scattered electron energy too, for more energetic particles originate 
more Cerenkov light, and so one can point back to their directions with greater accuracy. The probability to 
detect an electron at an angle a from the scattering direction is assumed to be a gaussian function 


dP{a, 




2 

da 


(2.9) 


where ctq = and is the empirically determined angular resolution of the detector for an electron 

with energy Cg- 

In Figure 3 one finds the the probability to detect an electron at an angle a from the original trajectory, as 
given by (2.9), for tg = 5, 15 MeV, and a plot of the angular resolution dependence on the electron energy for 
the Kamiokande-II detector a similar performance is expected for SK. In the case of a 5 MeV electron one 
has ~ 38°. 

Equation (2.9) tells us that one has a non-vanishing probability to detect a scattered electron just backward 
the supernova direction! This leads us to formulate a more reasonable definition for the forward observation cone 
than “the cone containing all v + e scattering events”. 

We verified that the model independent definition 

'dfwi^min) = \J '^'inaxi^rnin) + ('d^)(emm) (2-10) 

which can be read as the combination of the kinematic and detector spreads just like two independent errors, 
gives dfu, values such that no more than one v -\- e event is missed. For emin = 5 MeV one has dfw — 45°. The 
corresponding solid angle, II/u, = 27r [1 — cos{'dfw)], for emin = 5 MeV includes a 95% of the Vr events, but only 
a 15% of the total Dg signal. 

In Figure 4 one has the angular distribution of the tau-neutrino events, with respect to the supernova 
direction, considering both kinematic and detector spreads, while in Figure 5 we plotted the forward observation 
angle '&fw{f-min) as defined in (2.10). 
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Figure 1: Sketch of the directional analysis method. Electron (anti-)neutrinos coming from the supernova originate 
in the detector an isotropic signal (on the right). This one is used to infer the directional distribution for massless 
Vr (dashed line), which is to he compared with the observed signal (solid line), see bottom graph. 
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Figure 2: The angular distribution of Vr scattered electrons for two values of Cmin, and the maximum scattering 
angle iHmax os a function of the lower cut. Calculations refer to Burrows’s reference model, see Section 3.1 for 
details. 



a [deg] 


Figure 3: The probability to detect an electron at an angle a from the diffusion trajectory, for two values of Ce, 
and the experimental angular resolution in Kamiokande-II, from M- 
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Figure 4: The angular distribution of the Vt + e signal with respect to the supernova direction, including both 
kinematic and detector spreads, for two significative values of Cmin ■ Calculations refer to Burrows’s reference 
model (see Section 3.1). 



Figure 5: The forward observation angle as a function of the lower energy cut. 
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2.3 Supernova. 


Next, we have to describe how the Vi signal is derived from the Vf. one. For this purpose, we use three assumptions, 
supported by the results of several numerical simulations, (see ||l^) which connect the production features of 
electron and other flavor neutrinos in the cooling phase: 


1. The neutrino energy distributions, at any time and for each v species, are (approximately) thermal, thus 

specified by the temperatures of the neutrino-spheres and luminosities Lv^{t), 

The emitted neutrino rate spectrum dP/{de^ dt) at time t and energy Ci, is therefore: 

^ 0 _ ^I'x (P) _ ^ _ ('9 111 

de^dt “ l + ' 

where iFs = dxx^{l + e^)~^ ~ 5.68 

2. Energy equipartition among the six neutrino species holds at every time: 

L,At) = LuM = L.pt) ( 2 . 12 ) 


3. The temperatures of the six neutrino species are simply proportional at any time: 

(^) = (t) 

Tui {t) = a^. ■ (t) 

Arguments supporting the second assumption can be found in |l^ . Typical values of and are: 

= 0.5 - 0.9 

Op, = 1.1 - 1.9 

In what follows, we fix the ratio between and j?e-sphere temperatures at the central value Op^ 
we’ll keep Op. as a free parameter, in the range given by (2.14). 


(2.14) 
= 0.7, while 


2.4 Statistics. 

Finally, we introduce our criterion to estimate the evidence level of a non-vanishing mass. 

We define A^m([ti, ^ 2 ], [cmm, Cmax]) the expected number of events in the forward observation cone, in the time 
interval [^ 1 ,^ 2 ], for an energy window [cmm, f-max], if the Vt mass is m. 

Let us assume that tau-neutrinos have mass to 0: we expect the experiment to detect Nm events. The 
hypothesis that Vr have a vanishing mass and that the difference \Nra — AqI is due to a statistical fluctuation of 
the massless case can be tested by evaluating the quantity 

s(to, [ti,t2], [em*n,emax]) = (2.15) 

V Ao 

Clearly, one can vary [^ 1 ,^ 2 ] and [emin, £max] in order to optimize the analysis. The best indicator of a 
non-vanishing mass is therefore the quantity 

S'(to) = max {s(to, [tl,t2], [emm,emax])} (2-16) 

[^1 ,^2] j [^Tnin:^maxr 

where the time intervals are chosen arbitrarily within the observation window (first 100 seconds after bounce), 
with the only request that at least an event should be present in [^ 1 ,^ 2 ]- As an example, if 5'(100 eV) = 3 one 
has a 3cr evidence for a mass to = 100 eV, i.e. the Confidence Level is 99.7%. 

For brevity, we shall define as “detectable” those masses such that S{m) > 3, and we denote the minimal 
detectable mass as = min{TO : S{m) > 3}. 






Note that we are using only event counts, disregarding the possible additional information of spectral 
deformations. 

A sample of this statistical analysis is reported in Figure 6, where we plot the expected signal in the 
forward observation cone for m = 0 and m = 200 eV (for details on the model used for the calculation, 
see Section 3). For example, in the time interval [^ 1 ,^ 2 ] = [25,35] s, and for [emin,imax] = [15,25] MeV, one 
expects 1.3 events for the massless case, and 4.2 events for massive neutrinos; the evidence index is therefore 
s(200eF, [25, 35]s, [15, 25]MeF) = 2.5. We conclude that, in the time and energy intervals under examination, 
m = 200 eV is not sufficient to produce a clear effect. 


3 Results from theoretical models of supernovae. 

Till February 1987, numerical simulations were the only way to explore supernova neutrinos features. Most of 
the analyses just concerned the very initial phase of the emission {t ^ Is), mainly for implications on explosion 
mechanism, while only a few dealed with the long term cooling stage, which is relevant for us. For a review of 
simulations see for instance HI- 

In this section we use the simulations due to Burrows which describe the cooling phase, to evaluate 
the minimal mass m whose effects could be observed with the analysis method presented in Section 2. We also 
investigate the dependence of the results on the free parameters of the analysis. 

3.1 Models and general results. 

We performed a complete set of calculations using several supernova models taken from Burrows. In ref. [^ 
seventeen models are presented: disregarding “exotics” and black holes, the remaining ten models look consistent 
with data from S'iV1987A. 

Burrows’s data, that cover the first 20 seconds for 9^ luminosity and temperature, were extrapolated to later 
times by use of power-laws: 

A{t) = Ao(l + — 
ta 

where A denotes here both Ljj^ and and the parameters Aq, ta and ua are obtained from fits to late times 
behaviours (t > 10 s). We also introduced a time parameter ttr beyond which the neutrino luminosities vanish: 
this happens when the proto-neutron star becomes transparent to neutrinos. Typical values for ttr are (40 - 50) s 
[]l^ . Finally, for the other neutrino species, luminosities and temperatures are determined through relations (2.12- 
13). Unless otherwise stated, all results reported in the text, tables and figures, refer to the Super-Kamiokande 
detector for the following “default” set of parameters: 

= 1.5 

ttr = 45 s (3.1) 

D =10 kpc 

For each model, the main characteristics are listed in Table 2, together with the results: the optimal choices 
for time and energy windows, and the value of the minimal detectable mass, as defined in Section 2.4. The Vg 
event rate and the sensitivity to non-vanishing neutrino mass are shown in Figures 7 and 8 for a few representative 
models: the Most Favourable case (# 52), the Least Favourable case (# 62) and the FFference model (# 55), 
which very well accounts for S'V1987A data. 

We notice that the minimal detectable mass lies anyhow outside the cosmologically interesting region: even 
in the most favourable model we have msa — 120 eV. For the reference model one has m 3 a — 140 eV. 

When looking at the dependence on the model characteristics, we find that in general: i) for a given equation 
of state (EOS), mass effects are least pronounced for more massive cores; ii) for a given core mass, neutrino mass 
effects are weaker for the soft EOS. 
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Those features can be understood, at least qualitatively. We remark that, for a given EOS, the star radius 
decreases with mass, so that in more massive cores matter is more compressed. Also, for a given mass, softer 
EOS allow higher compression. Thus higher masses and/or softer EOS imply higher densities, and consequentely 
higher temperatures and opacities. More energetic neutrinos are produced, so that i/'T--mass induced delays are 
shorter. Furthermore, and more important, the longer cooling phase results in a higher background Vf. flux at 
later times. The increase of the cooling time due to mass and EOS is shown in Figure 7: ME and RE models 
have the same EOS but different masses, whereas LF and RE models have the same mass but different EOS. 

3.2 Dependence on the model parameters. 

In order to investigate the dependence of the previous results on the emission parameters, we repeated our 
calculations by varying one parameter at a time, while keeping all the others constant at the values of (3.1). The 
results quoted in the following are for the reference model, but we obtained similar behaviours for the others. 

We first varied the parameter defined in Section 2.3 as the ratio of Vi- to -sphere temperatures, 
repeating the analysis for hve different values in the range given by (2.14). 

The results of this search are summarized in Table 3, where we report, for each value of the parameter, the 
minimal detectable mass, together with the number of v-r events and the characteristic delay time for m = 150 eV. 
In Figure 9 we plot the evidence index S{m) for the central value and the two extrema of the ai,^ range. The 
minimal detectable mass depends rather weakly on the z/^-sphere temperature. The point is that when varying 

two competing effects arise. For example, when an^ is increased: i) the z^T-sphere temperature is higher, 
leading to more energetic neutrinos, and so to shorter delays; ii) the number of events increases, leading to 
a higher statistics. Our calculations essentially show that these effects almost balance each other, so that small 
variations of a^i do not affect substantially the sensitivity to a non-vanishing mass. 

Concerning the dependence on the transparency time, we performed three calculations of the minimal 
detectable mass for ttr around the value given by [^. Again the results are essentially stable, see Figure 
10, where we plotted the evidence index S(rn) for the inquired values. Clearly, the more sudden the transparency 
occurs, the more evident the delayed neutrinos will be: should we take the (maybe) unphysical choice ttr = 30 s, 
we find m^rr — HO eV. 

Finally,.... 

3.3 Best choices for the analysis parameters. 

These are particularly important dependences to investigate, because we deal now with detection parameters, 
which we can change in order to optimize the analysis. 

It was already noted 0 that an upper energy cut enhances the mass sensitivity of the analysis method. This 
is due to the fact that the isotropic signal occurs at a higher energy with respect to the directional one. In fact in 
the dominant isotropic process, i.e. De capture on a proton, the emitted positron takes almost all of the neutrino 
energy, while in the Vr process, namely scattering off an electron, the average energy of the outgoing particle is 
roughly ^ ^e,y. We stress that by means of e^ax we cut high energy i^r events, which have short mass-induced 
delays, so that the loss in sensitivity is low. We should be aware, anyway, that this dynamical effect is partially 
balanced by the higher Vr energy at the emission. 

Concerning the lower energy cut that we introduced, many competing effects are present, leading to a more 
complex situation. Basically, at high values of Cmin the number of Vt events get small and, most of all, we deal 
only with high energy neutrinos, for which the flight-time delay is shorter. On the other hand, as we said in 
Section 2.2, the forward observation angle is a decreasing function of the lower cut (see eq. 2.10). Rising emin 
leads therefore to a narrower cone, in which the ratio between directional and isotropic Ve signal becomes more 
favourable. 
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Figure 6: An example of the statistical analysis: the dashed line refers to the expected distribution in the forward 
observation cone for the case m = 0, while the solid one is for m = 200 eV, see text for explanations. Burrows’s 
“reference” model described in Section 3 was used for this example. 


Model 

EOS 

M, [Mq] 

Mf [Mo] 

Etot [foe] 


[tl,t 2 ] [s] 

mscr [eV\ 

52 

stiff 

1.2 

1.2 

152 

[10,20] 

[25,85] 

120 

53 

stiff 

1.3 

1.3 

171 

[10,20] 

[30,95] 

130 

54 

stiff 

1.4 

1.4 

191 

[10,20] 

[30,100] 

130 

55 

stiff 

1.3 

1.5 

228 

[ 10 , 20 ] 

[40,95] 

140 

56 

stiff 

1.3 

1.6 

262 

[10,20] 

[45,100] 

130 

57 

stiff 

1.3 

1.8 

326 

[10,20] 

[40,95] 

130 

59 

soft 

1.2 

1.2 

158 

[10,20] 

[35,100] 

140 

60 

soft 

1.3 

1.3 

178 

[10,20] 

[40,100] 

140 

61 

soft 

1.3 

1.4 

205 

[10,20] 

[40,100] 

140 

62 

soft 

1.3 

1.5 

235 

[10,20] 

[40,100] 

140 


Table 2: Characteristics and results for Burrows’s supernova models: the columns indicate the equation of state 
(EOS), the progenitor core mass Mi, the proto-neutron star mass Mf (after accretion), the total emitted energy 
Etot (1 foe = 10 ®^ erg), the optimal energy [emim^max] and time windows in order to probe mass effects, 

and the minimal detectable mass m^a- 
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Figure 7: The isotropic Ve signal in the whole SK detector, for the least/most favourable model (LF/MF), and 
for the “reference” model (RE), see text. 
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Figure 8: Mass sensitivity for a few Burrows’s models: same notations as in Figure 7. 
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Figure 9: Mass sensitivity for several values of the Vi-sphere temperature, in Burrows’s reference model. 




At{ei, = a,,^■15MeV) [s] 

[^min-; ^max] [AfeV^] 

[tut2] [s] 

msa [eV] 

1.0 

4.93 

52 

[10,20] 

[35,95] 

140 

1.3 

11.8 

31 

[10,20] 

[35,100] 

140 

1.5 

16.4 

23 

[10,20] 

[40,100] 

140 

1.7 

16.8 

18 

[10,20] 

[40,100] 

140 

2.0 

17.5 

13 

[10,20] 

[40,100] 

150 


Table 3: Mass sensitivity dependence on the Vi-sphere temperature in Burrows’s reference model, with ttr = 45 s. 

is the number of Vr events in the forward cone and At is the characteristic delay time for a neutrino in the 
explorable mass range (m = 150 eV) and for D = 10 kpc. 


ttr [ 5 ] 

^u-T-ittr ^ tuiax) 


[^ 1 ,^ 2 ] [s] 

mso- [eV] 

30 

6.2 

[10,20] 

[30,70] 

110 

40 

3.9 

[10,20] 

[35,95] 

130 

45 

3.2 

[10,20] 

[40,100] 

140 

50 

2.5 

[10,20] 

[40,100] 

140 


Table 4: Mass sensitivity dependence on the transparency time, in Burrows’s reference model, with a^i = 1-5. 
Alv^{ttr tmax) is the number of Vr events expected to occur after the transparency, in the case of m = 150 eV 
and for D = 10 kpc. 
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Figure 10: Mass sensitivity for several values of the transparency time, in Burrows’s reference model. 



Figure 11: Mass sensitivity for several values of the supernova distance, in Burrows’s reference model. 
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In Table 5 we report the full results of our analysis on the reference model for six different energy windows, 
while in Figure 11 we plot the mass sensitivity for three significative choices. 

As we can see from Figure 11, the best choice depends on the explored mass region, because of the m 
dependence in the delay formula (1.1): the larger is m, the higher should be the lower cut. For our explorable 
range m = 100 — 200 eV the choice emin = 10 MeV ; e^ax = 20 MeV gives the best sensitivity for all models and 
for all values of the emission parameters ai,^ and ttr- 

With regard to the time window, as we can see from Table 2 the best choice again depends on the explored 
mass: larger m mean more delayed neutrinos, which are to be observed at later times. For our explorable range 
the choice ti = 40 s; t 2 = 100 s seems in general the most appropriate. 

3.4 This section’s conclusions. 

We conclude, for the case of a Burrows-like supernova near the galactic center, that: 

1. there is no way of getting evidence of m ^ 120 eF; 

2. this conclusion is essentially independent of the value of the uncertain emission parameters and ttr', 

3. for the explorable mass range, the optimal energy window is Ce ~ [10, 20] MeV. 

4. for the explorable mass range, the optimal time window is t ^ [40,100] s. 

4 Results from phenomenological models of supernovae. 

So far we considered theoretical models for neutrino emissions, essentially based on numerical simulations. In 
order to investigate further the potentialities of the analysis method, in this section we consider phenomenological 
models, built so as to reproduce the data from 5^19874. 

The mere 19 events of the 1987 supernova were not sufficient to allow a detailed reconstruction of the emission 
features, and several models were obtained by fitting the data with different analytical functions. We selected 
two of these models, belonging to significative classes, to explore a wider range of possibilities. 


4.1 Power-law cooling. 

In this framework, we first consider a power-law parametrization, defined by the relations: 

r T{t) = TO (1 + ^)“” 

\ m = L0{i + i)-^- 

Note that, being L{t) oc in this model the radii of the neutrinospheres are assumed as constant. 

We use as best fit parameters for SN1987A those obtained by Bludman & Schinder |l^: 


'T'O 

~ 4.20 

MeV 

Ll 

~ 19.8 

foe/s 

T 

~ 2.78 

s 

n 

~ 0.4 



(4.1) 


(4.2) 


Again we used relations (2.12-13) to derive luminosities and temperatures for the other neutrinos species. 

As we can see from Figure 12, where we plotted the Dg event rates in SK for (4.1) and for Burrows’s reference 
model, the behaviour of the power-law cooling is similar to that of theoretical models of Section 3. 

With the values of (4.2) and for the default set of parameters (3.1), we obtained the minimal detectable mass 
771-30- — 110 eV. This result is consistent with what found for Burrows theoretical models, and is again outside 
the cosmologically interesting range. 
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D [kpc\ 


At(e^ = 20MeF) [s] 


[h,t 2 ] [s] 

m^a [eV\ 

8 

25.5 

14.4 

[10,20] 

[40,100] 

140 

10 

16.4 

18.0 

[10,20] 

[40,100] 

140 

12 

11.4 

21.6 

[10,20] 

[35,100] 

130 

50 

1.3 

90.0 

[5,100] 

- 

- 


Table 5: Mass sensitivity dependence on the supernova distance, in Burrows’s reference model, with = 1-5 and 
ttr = 45 s. is the number of Vr events and At is the typical neutrino time-delay in the case of m = 150 el^. 


6 

> 

> 

5 

4 

^ 3 

ui 

2 

1 

0 

0 50 100 150 200 

m [eV] 

Figure 12: Mass sensitivity for significative choices of the energy window [emin,emax\, in Burrows’s reference 
model. 
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N,.jNtot [%] 

'dfw 

N.jNtot [%] 


[^ 1 ,^ 2 ] 


[MeV] 

(all SK) 

[deg] 

(in 

= — ^min) 

[s] 

[ell 

[7.5,15] 

0.91 

37 

6.4 

210 

[40,100] 

140 

[10,20] 

0.66 

32 

6.5 

120 

[40,100] 

140 

[10,100] 

0.52 


5.4 


[30,100] 

190 

[15,25] 

0.48 

27 

7.0 

52 

[20,80] 

160 

[15,100] 

0.43 


6.4 


[20,90] 

170 

[20,100] 

0.39 

24 

6.8 

29 

[20,90] 

180 


Table 6: Mass sensitivity for significative energy windows, in Burrows’s reference model. Calculations are for 
ttr = 45 s and ai,^ = 1.5. Atmax is the maximum delay for a neutrino in the explorable mass range (m = 150 eF ), 
corresponding to the case of minimal detected energy, for D = 10 kpc. 
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Figure 13: Total signals in the whole detector for exponential cooling (dotted), power-law cooling (dashed), and 
Burrows’s reference model (solid). 
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4.2 Exponential cooling. 


A simpler analytic parametrization is obtained assuming that the neutrinosphere temperatures fall exponentially, 
while the radii are again kept constant. One has this way 


T{t) = 

L{t) = LOg-t/r 


The best fit parameters for SN1987A are taken from Loredo & Lamb 




(4.3) 


'T'O 

~ 4.16 

MeV 


Ll 

~ 12.9 

foe/s 

(4.4) 

T 

~ 4.6 

s 



With these choices one obtains a minimal detectable mass mso- — 70 eV. The exponential cooling, although 
not clearly distinguishable from the power-law on the grounds of the SA^IOSTA data, predicts at later times 
{t ^ 20 s) a substantially reduced Pe hux (see again Figure 12). This behaviour can account for the sensitivity 
to smaller masses: in fact, with a lower 9^ background, the presence of delayed massive neutrinos is much more 
evident. 

Actually, there is no physical reason for assuming an exponential law, and theoretical models generally 
predict a slower decrease of neutrino luminosity at late times Furthermore, when compared with S'iV1987A 
data, exponential cooling models yield fits which are worse than those of power-law models . In conclusion, it 
seems to us that exponential cooling models are not very reliable for predicting late times neutrino luminosities. 


5 Comparison with Krauss’s results. 

Krauss et al. [^, once developed their analysis method, also evaluated the minimal detectable mass. They 
found that a mass as small as 50 eV is detectable at 99% C.L. with Super-Kamiokande, for a medium luminosity 
burst from a SN at D = 10 kpc. 

This result is more optimistic than we found by using Burrows’s theoretical models and/or 5'iV1987A 
phenomenological models (we recall, e.g., that for Burrows’s most favourable model one has toso- = 120 eV at 
99.7% C.L.). 

This discrepancy can’t be ascribed to the differences in the analysis procedures. We have verified that, with 
our approach, one obtains for the middle model of (n. 17) a minimal detectable mass mso- — 48 eV. 

Actually, in the supernova models of after an initial accretion phase {t ~ 0.5 s), the luminosities of all 
neutrino flavors are assumed to decay exponentially with a very fast rate {jl ~ 1 s). This results in a quite short 
tail of the event distribution, and so the signal of delayed Vt is much more evident, as we already noted for the 
5^^19874 exponential model. 

This feature is evident if we compare a typical model from [jl^ with one of those we used in the foregoing. In 
Figure 13, for instance, we plotted the isotropic (scaled to ’df^) and directional event rates both in Krauss’s model 
17 and in Burrows’s model 54, which give fairly the same event number. As we can see, with just m = 100 eV, in 
Krauss’s model the directional signal outruns the isotropic “background” after a few seconds, while in Burrows’s 
model this occurs only at a very later times (t > 30 s). 

It looks to us that the fast luminosity decay assumed in jl^ is not justified since: i) phenomenological 
models with exponential cooling yield longer decay times when adapted to fit 5'A^19874 data, see §M ; ii) also 
theoretical models predict longer decay times, see and the previous discussion. 

^We remind that cooling of a degenerate gas by black-body radiation yields a power-law decrease of temperature: 
^(^)=T“(l-b|)■^^^ 
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6 Upper bounds. 


The situation may look somehow discouraging; nevertheless we should remember that the present experimental 
limit to Vt mass is tens MeV. Therefore even in the most probable case of non-observation of mass effects, we 
will however significantly lower this limit. 

In order to obtain actual upper bounds, anyway, we have to introduce a more stringent statistical criterion 
than the one we used to define our “detectable” mass. With our definition of Section 2.4, in fact, we calculated 
the minimal value of m expected to produce a 3cr effect, checking against statistical fluctuations of the massless 
case. We should now take into account fluctuations of both massive and massless case. 

Therefore, the new best indicator is the quantity 

o! \ — \^m ~ fVol /r T\ 

b{m) = max ^ (b- 1 ) 

^TYl + No 

As an example, if 5(100 el^) = 3 we can claim the upper bound m < 100 eV at a 3cr level, i.e. with a 99.7% 
Confidence Level. 

In practice, Nm is generally (much) larger than Nq, so that (6.1) can be replaced to a good approximation 
with: 

S{m) = max l-^m -^ o[ 1-0 2 ) 

kl.t2].k™m.emax] ^/Wrn 

This can be read very much like (2.15), i.e. assuming that neutrinos are massless, the observer will detect Nq 
events (or a number close to it). Eq. (6.2) tells then how many sigma’s the result is out of the expectation for 
massive neutrinos. 


Model 

TOScr [eV\ 

[^mint ^max] 

[U,t 2 ] [s] 

mso- [eV\ 

Burrows RE 

140 

[15,100] 

[15,100] 

200 

Burrows ME 

120 

[ 10 , 20 ] 

[15,95] 

140 

Power-law 

no 

[ 10 , 20 ] 

[10,85] 

140 

Exponential 

70 

[ 10 , 20 ] 

[10,60] 

90 

Krauss 17 

50 

[ 10 , 20 ] 

[5,15] 

70 


Table 7: The minimal excludable mass for a few significative models of supernova. 


The results for two significative theoretical models (Burrow’s reference and most favourable models - see 
Section 3) and two phenomenological models (exponential and power-law cooling - see Section 4) are listed in Table 
6 . Using the same notation of Section 2.3, we define the minimal excludable mass as ffisa = min{TO : S{m) > 3}. 
Clearly a larger mass is required to probe an upper bound rather than a simply detectable effect, for one has 
Nm > Nq ^ S{m) < S{m). Anyway, one can see that even in the case of non-observation of mass effects, we 
will able to lower the present upper bound on Vr mass by 5 orders of magnitude! 

One should be aware, anyway, that the bounds of Table 6 (and also the previous results on the minimal 
detectable mass) may not hold if the mass is too large. In this case, in fact, one has a too broad event 
distribution, see for istance [ po| . A very crude estimate of the maximum mass mmax explorable by delay-based 
analyses can be obtained by requiring that the rate of delayed Ur events keeps greater than the background 
rate. For SK one finds roughly rumax ^ 0(10 keV). Our results should be read therefore as lower extrema of 
an excluded mass region: in the case of non-observation of mass effects one could exclude m in the window 
200 eU-20 keV. 
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7 Conclusions. 


Let us summarize the main points of this paper: 

i) We performed an extensive investigation of the sensitivity to non-vanishing Vr mass in a large water 
Cerenkov detector, developing the analysis method introduced by Krauss et al. jf^. As the most important 
point is the dependence of the results on the supernova model, we analysed several theoretical models from 
numerical simulations of Burrows [|^ and phenomenological models based on SN1987A data @1- 

ii) We determined optimal values of the analysis parameters (energy cuts and time window) so as to reach the 
highest sensitivity to a non-vanishing Vr mass. 

iii) Our conclusion is that the minimal detectable mass is generally just above the cosmologically interesting 
range, m ^ 100 eV, see Table 2 for details. 

iv) Differences with respect to 0, which claimed sensitivity to signihcantly smaller masses, are ascribed to 
the shorter decay time of neutrino luminosities assumed in pT[ |, contrary to theoretical models as well to 
SN1987A phenomenological models. 

v) For the case that no positive signal is obtained, observation of a neutrino burst with SK will anyhow lower 
the present upper bound on Vt mass (m ^ 24 MeV Q) to few hundred eV, see Table 6. 

A Appendix 

A.l Super-Kamiokande characteristics. 

In Table 7 we summarize the main characteristics of the Super-Kamiokande detector (the fiducial mass, the 

threshold energy and detection efficiency, and the background rate spectrum) as we used in our analysis. 



^ We refer to the hducial mass for supernova events. 


Table 8: Super Kamiokande characteristics. 


A.2 Cross sections. 


1. capture on a proton: 


da. 




dCe 






(1 -f 3a^) (1 -l- Swm) ( 1 





Q) (A.l) 


where a ~ —1.26 is the axial-vector coupling constant, Q ~ 1.293 MeF is the neutron-proton mass 
difference, S^m — —3.3 • 10“^(e;y — Q 12)1 MeV is the weak-magnetism correction, and 


4 1 


(To = - ■ 
TT 


ch)' 


-Gl 


^)2 ~ 1.76 • lO"^"* cm? 
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2. scattering off an electron: 


duu 


1 (To 


dCe 2 nrieC^ 

where the constants A and B have values: 


A^ — 




(i -b sin^ 6wY 

~ 0.536 

Ve 

sin^ 9w 

~ 0.0538 , 

Ve 

(-i + sin^ 9w)'^ 

~ 0.0719 , 

Vi 

sin"^ 9w 

~ 0.0538 , 

Vi 

sin"' 9\y 

~ 0.0538 , 

Ve 

(i -b sin^ 9w)^ 

~ 0.536 , 

Ve 

sin"' 9w 

~ 0.0538 , 

Vi 

(-i + sin^ 9w)‘^ 

~ 0.0719 , 

Vi 


Brr = 


3. capture of a v^. on an oxygen nucleus: 

The cross section of this process can be roughly approximated to plf 


da. 


L'eO 


~ 0.16(To 


where 


dce \ irieC 

13 MeV is an effective reaction threshold. 


9{ei, - e) • i5(ey - e - Ce + mff 


4. capture of a Ve on an oxygen nucleus: 
We use the approximation given by [Dl 


da. 


UeO 


~ 0.074(To 


dte \ rneC 

where the reaction threshold is again e ~ 13 MeV. 


6{ey - e) • 5[ev - e - Ce + me) 


(A.2) 


(A.3) 


(A.4) 


A.3 Prom emission to detection. 

Being df’N^^/dcv dt' the rate spectrum of Vx emitted at time t' with energy as defined in Section 2.3, the 
corresponding differential flux at the detector at time t is 

d‘^4>^ffe^,t) 1 f d‘^N^ffe^,t') 

d,„dt =^J'“ d,„dt' (A,5) 

For each detection process one obtains therefore the event rate dNb/dt as 

da^b{eu,ee) 
lib 


dNb{t) r— ( ^ r n 


' dt ^ ^ ^ 

where da^b/dcg is the differential cross section of the reaction under examination, Ub is the target number, and 
ri{ee) is the electron/positron detection efficiency. 

For the natural background one has instead the time-independent rate 


B{e 


min 1 ^max ) 


dce b{ee)r]{ee) 


(A.7) 


where &(ee) is the empirically determinated background rate spectrum. 
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Figure 14: Isotropic (scaled to Dfw) and directional signals in Krauss’s model n. 17 and in Burrows’s model n. 
54 (same event number); a m = 100 eV was assumed. 
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